Hydrogenated Graphene Nanoribbons for Spintronics 
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We show how hydrogenation of graphene nanoribbons at small concentrations can open new 
venues towards carbon-based spintronics applications regardless of any especific edge termination 
or passivation of the nanoribbons. Density functional theory calculations show that an adsorbed H 
atom induces a spin density on the surrounding n orbitals whose symmetry and degree of localization 
depends on the distance to the edges of the nanoribbon. As expected for graphene-based systems, 
these induced magnetic moments interact ferromagnetically or antiferromagnetically depending on 
the relative adsorption graphene sublattice, but the magnitude of the interactions are found to 
strongly vary with the position of the H atoms relative to the edges. We also calculate, with the 
help of the Hubbard model, the transport properties of hydrogenated armchair semiconducting 
graphene nanoribbons in the diluted regime and show how the exchange coupling between H atoms 
can be exploited in the design of novel magnetoresistive devices. 



I. INTRODUCTION 

There is a widespread consensus on the large pote- 
tial of graphene for electronic applicationsiiS. Theoret- 
ically, graphene also holds promise for a vast range of 
applications in spintronics, although clear evidence of 
magnetic graphene is, however, elusive to date. In a 
broad sense, two factors may account for this elusive- 
ness. First, the fact that hydrocarbons of high spin are 
known to be highly reactive and, unless fabricated or syn- 
thesized under very clean and controled conditions, they 
will likely bind surrounding species with the concomit- 
tant disappearance of magnetisnt^. A second reason re- 
lates to the fact that the ground state of graphene is near 
an interaction-driven phase transition into an insulating 
antiferromagnetic stated. This underlying antiferromag- 
netic correlations prevent the magnetic moments, even if 
they develop, from ordering ferromagnetically and pre- 
clude the possibility of observing hysteresis in standard 
magnetic measurements. Notwithstanding, a few reports 
of magnetic graphite^ and graphene^ can be found in re- 
cent literature. 

Most of recent theoretical ideas for graphene based 
spintronics applications are rooted on the magnetic prop- 
erties of nanoribbons!"— or nanographenesi^"— with 
zigzag edges. All these proposals assume a very par- 
ticular edge hydrogenation where H atoms passivate 
the a dangling bonds, leaving all the tt orbitals un- 
saturated and carrying the magnetic moments. How- 
ever, this is just one out of many possible edge realiza- 
tions which range from H-free self-passivationi^ to full 
H passivationii. According to the work of Wassmann et 
alJi, relatively low H concentations at room tempera- 
ture suffice to completely passivate the edges, including 
the edge tt orbitals responsible for the magnetic order. 
The self-passivated or reconstructed zigzag edges^^ are, 
in fact, among the least energetically favorable of all, 
although, interestingly, have been recently observed by 



transmission electron microscopylS. 

In the light of the present controversy on the actual 
possibilities of ever encounter zigzag magnetic edges, we 
propose in this work an alternative to edge-related spin- 
tronics in which to exploit the recently shown controled 
hydrogenation of graphene-^. The key factor here is that 
adsorption of atomic H in the bulk of graphene is ac- 
companied by the appearance of a magnetic moment 
of IfiB localized on the tt orbitals surrounding each H 
atoni22,. These magnetic moments interact with one an- 
other ferromagnetically or antiferromagnetically, depend- 
ing on whether their respective adsorption sublattices 
(usually labeled A and B) are the same or not3"— . Sta- 
tistically speaking, a sublattice compensated H coverage 
is expected unless adsorption on one sublattice is privi- 
leged, e.g., by the substrate. To date, however, there is 
no evidence that such an uncompensated coverage can be 
achieved. In the more likely compensated case an overall 
antiferromagnetic alignment with a total spin 5' = is 
thus energetically favored over a ferromagnetic one with 
5 > 0. As a proof of principle and since we are inter- 
ested in the diluted regime, we consider in this work the 
fundamental problem of two H atoms. These are cova- 
lently bonded to the surface of a semiconducting arm- 
chair graphene nanoribbon (AGNR) on different sublat- 
tices [see Fig. (Ha)]. Here the a bonds of the edges are 
fully passivated with H so that they arc irrelevant at low 
energies. 

As shown in Sec. IIIIl after a brief introduction to 
the theoretical basics presented in Sec. [TTl the magnetic- 
field driven ferromagnetic (F) state, where the H-induced 
magnetic moments are aligned by the field [as shown in 
Fig. IHb)], can present a different resistance from that 
of the natural antiferromagnetic (AF) zero-field state [as 
that in Fig. [H^a)]. Two different cases are discussed: Infi- 
nite semiconducting AGNR's and finite ones connected to 
conductive graphene [see schematic piture in Fig. (Hd)]. 
The differences in conductance between the F and the AF 
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FIG. 1: (a) Armchair graphene nanoribbon of width A'' (where 
A'' is the number of dimer hnes) with two H atoms (shown in 
green) adsorbed in the middle of the ribbon at a distance d 
from each other. The inset shows a detail of the adsorption ge- 
ometry, (b) Pictorical representation of the one-orbital tight- 
binding model where the presence of H atoms is simulated by 
removing sites in a head-to-head configuration, (c) Same as 
in (b), but for the opposite ordering (tail-to-tail), (d) Same 
as in (b), but for a finite semiconducting AGNR connected to 
metallic nanoribbons at the edges. 




(b) 

FIG. 2: Pictorical view of the antiferromagnetic state (a) and 
the magnetic-field driven ferromagnetic state (b) where the 
magnetic moments localized around two hydrogen atoms are 
depicted by red arrows (the orientation of the arrows with 
respect to the graphene plane is arbitrary since spin-orbit 
coupling is neglected here). 



states can be substantial and translate into a magnetore- 
sistive response as large as 100% for distances between 
the H atoms of the order of few nanometers. Practical 
implications of these results are discussed in Sec. ITVland 
a brief summary is presented in Sec. |V] 



II. COMPUTATIONAL DETAILS 

For the calculation of the electronic structure of hydro- 
genated AGNR's we use both ab initio techniques within 
the local spin density approximation (LSDA), aided by 
the CRYSTAL03 package^^., and a one-orbital (tt) first- 
neighbor tight-binding model where the electronic repul- 
sion is treated by means of a Hubbard-like interaction U 
in the mean field approximation: 

i,j i i 

. (1) 

The first term represents the kinetic energy with first 
neighbors hopping t between tt orbitals. The remain- 
ing terms account for the electronic interactions where 



•hijj — cf^Ciri are the number operators associated to 
each TT orbital with spin 77. These two different lev- 
els of approximation to the electronic structure have 
been shown to yield similar results for nanoribbons^^ and 
nanographeneai^ where a full passivation by H of the <j 
bonds is assumed. 

The bulk adsorption geometry of a H atom has been 
obtained by relaxing the C atom bonded to the H and the 
nearest C atoms until the characteristic sp^ hybridization 
is obtained [see inset in Fig. (Ua)]. The bonding between 
a H atom and a C atom results in the effective removal 
of the TT orbital from the low energy sector, so that the 
H adsorption is simulated by simply removing a site in 
the one-orbital mean-field Hubbard model [see Fig. [Ifb- 
d)]. Our results, as shown below, provide further support 
to the use of the Hubbard model in graphene systems, 
as an alternative to the computationally more demand- 
ing LSDA and extend the range of applicability of Lieb's 
theorem^^ to a wider set of situations. 

To calculate the transport properties we use the 
standard Green's function partitioning method as im- 
plemented, e.g., in the quantum transport package 
ALACANT^^. To this purpose, the infinite system is 
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divided into three parts, namely a central region (C), 
containing the H atoms, which is connected to the right 
(R) and left (L) semi-infinite clean leads. The Hamil- 
tonian matrix describing the whole system is then given 

by 



H = Hc + Hr + Hl + Vlc + Vrc 



(2) 



where He, and He. are the Hamiltonian matrices of 
the central region, the left and the right lead, respec- 
tively. Vlc and Vrc represent the coupling between 
the central region and the leads. In general, the non- 
orthogonality of the basis set must be taken into account 
when writing the Green's function of the central region: 



Gc{E) = [ESc -He- ^l{E) - ^r{E)]- 



(3) 



The self-energies of the left (Sl) and right (S_r,) leads 
account for the influence of these on the electronic struc- 
ture of the central part and Sc is the overlap matrix in 
this region. For the calculation of the conductance we 
use the Landauer formula, G{E) ~ ^T{E). The trans- 
mission function T can, in turn, be obtained from the 
expression; 



T{E) = ^Tr [gUe)Tr{E)Gc{E)Tl{E) 



(4) 



with T j- 



Notice that since we 



are interested in coUinear magnetic solutions, all the ma- 
trices carry the spin index rj, which we have not made 
explicit in previous equations. All the terms in Eq. [2] 
must be obtained self-consistently either from a periodic 
boundary condition calculation, e.g., using CRYSTAL03 
in the case of the LSDA calculations, or following the 
methodology in Ref.^ in the case of the Hubbard model. 
The Fermi energy is set to zero and, in both cases, global 
charge neutrality in all regions is imposed by shifting the 
onsite energies as necessary. 



III. RESULTS 

A. Energetics 

We first examine the energetics of the F {S = 1) and 
AF {S = 0) states as a function of the mutual distance 
d between H atoms. We first choose a semiconducting 
AGNR of width = 9, where N is the number of dimer 
lines across the ribbon, and restrict ourselves to the case 
of H atoms placed on different sublattices [see Fig. [H^a)]. 
The reason for this choice is three-fold. The AB (or BA) 
configurations are always energetically preferred to the 
AA or BB configurations for similar distances between H 
atoma^Si. Second, and most importantly for the purpose 
of this work, the magnetic state of the AB (or BA) config- 
uration can be tuned by a magnetic field. Furthermore, 
as briefly mentioned in the introduction, even if energetic 
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FIG. 3: (a) Total energy referred to Eq [lowest energy in case 
(b)] as a function of the distance between H atoms for the fer- 
romagnetic (dashed) and the antiferromagnetic (solid) state 
when placed in the middle of the ribbon. Upper inset: Pic- 
ture of a nanoribbon with two H atoms. Lower inset: Energy 
difference between both states and extrapolation to large dis- 
tances (solid line). The vertical line in the inset denotes the 
distance above which the energy difference becomes less than 
1 meV. (b) Same as in (a), but for both H atoms near the 
same edge, (c) Same as in (a), but for H atoms placed on op- 
posite edges. Here the width of the ribbon is larger (A^ = 13) 
and Eo is the minimum energy among all the AF solutions. 
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considerations are left aside, the AB (or BA) configura- 
tion represents the simplest case of a random ensemble 
of H atoms which, in average, will equally populate the 
two graphene sublattices. 

In Fig. ISlJa) we plot the LSDA total energy for H 
atoms in the middle of the AGNR as a function of d. 
The S = state is always the ground state. This implies 
antiferromagnetic coupling, except at short distances 
for which the local magnetization, quantified through 
^ = vanishes altoghether [see Fig. HJ^c)]. 

The quenching of the magnetization is easily understood 
in terms of the formation of a spin single1j21. At the 
minimum distance for which S 7^ 0, the energy differ- 
ence presents a maximum and decays exponentially for 
larger d [see inset in Fig. [31[a)]. As d 00 the spin 
clouds do not interact anymore and both F and AF so- 
lutions tend to have the same energy. In summary, for 
any distance between H atoms the ground state presents 
S = 0, following Lieb's theoreni'^^, but the overall spin 
texture strongly depends on their mutual distance. We 
note that it also depends on the ordering (AB or BA) 
of the H atoms. Whereas in bulk the spin cloud associ- 
ated to the H atom would be invariant under rotations 
of 120 degrees, in a nanoribbon there is a preferential 
direction along the ribbon axis which is different for H 
atoms located on A and B sublattices [see Fig. [IJb-c)]. 
We refer to the preferential direction as the head and the 
tail to the opposite one. Thus, a head-to-head coupling 
(AB) is expected to be much stronger than a tail-to-tail 
coupling (BA). The magnetization clouds are shown for 
an AB (or head-to-head) case in Fig. UJ^a-b). It is easy 
to appreciate the strong directionality just alluded to. 
As a consequence, when reversing the ordering of the 
H atoms to a BA (or tail-to-tail) configuration [see Fig. 
[IKc)], these do not couple magnetically at any distance, 
except when in very close proximity for which E = 0. All 
these results are similar to the ones reported using the 
one-orbital mean-field Hubbard model^l. 

In Fig. in^b) we present the LSDA total energy of the 
F and AF states when the H atoms are placed near one 
of the edges. Again the ground state is the AF state for 
any distance. The proximity to the edge decreases the 
localization length of the spin texture, increasing E [see 
Fig. Hfd-e)], and thereby, decreasing the critical distance 
below which the magnetization disappears [see Fig. IDJf)]. 
The perturbation of the edge modifies the spin texture 
around each atom and, contrary to the previous case, 
when reversing the ordering to a tail-to-tail configura- 
tion, the H atoms couple magnetically in a finite range 
of distances (see below). All energies, including those de- 
picted in Fig. [3lja), are referred to Eq, the lowest energy 
solution for H atoms close to the edge, corresponding 
to the smallest distance. Note that for the same distance 
between H atoms the total energy is always smaller when 
these are closer to the edge, which reflects that the bind- 
ing energy is higher there by approximately 1 eV— . 

Finally we present in Fig. ^c) the case where the 
H atoms are placed on opposite edges for an = 13 
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FIG. 4: Magnetic moments (whose magnitude is represented 
by the size of the circles) on individual C atoms when the 
H atoms are placed head-to-head in the middle of an A'^ = 9 
armchair graphene nanoribbon at a distance d = 32.0 A [(a) 
ferromagnetic state and (b) antierromagnetic state] and at a 
distance d = 15.4 A close to the edge [(d) ferromagnetic state 
and (e) antiferromagnetic state] . Panels (c) and (f ) show the 
magnetization (see definition of S in text) as a function of the 
distance between H atoms. 



semiconducting AGNR. Placing the atoms on different 
edges allows us now to explore the energetics of differ- 
ent magnetic coupling orientations (tail-to-tail for d > 
and head-to-head for d < 0) at a reasonable computa- 
tional cost (d is now the longitudinal distance between 
H atoms). As expected, due to the strong anisotropy 
of the spin texture, the magnetic coupling strongly de- 
pends on the ordering of the atoms and not only on their 
mutual distance. While for d < the energy difference 
between the AF and F states is large, for d > this dif- 
ference decreases substantially and does not depend too 
much on d for the values considered. This asymmetry 
is easily understood since for d > 0(d < 0) the spin tex- 
tures approximately couple in a tail-to-tail(head-to-head) 
manner. As shown below, this may have important ex- 
perimental consequences. 



B. Magnetoresistance 

We now turn our discussion to the implications these 
results may have on the electrical transport. Under the 
influence of a magnetic field, the hydrogenated AGNR 
behaves like a diluted paramagnetic semiconductor^^ for 
small concentrations of H. At large concentrations, when 
the spin density is zero everywhere, the infiuence of the 
field can only give rise to a minor diamagnetic response. 
At intermediate concentrations, where the magnetization 
clouds induced by the H atoms interact with each other, 
one can switch from the AF to the F state by apply- 
ing a sufficiently strong magnetic field. In analogy with 
the II2 molecule, where switching from the singlet to the 
triplet state modifies the orbital part of the wavefunc- 
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FIG. 5: Spin resolved transmission as a function of energy for 
the ferromagnetic state of an armchair graphene nanoribbon 
of width N = 9 with two H atoms adsorbed in the middle at 
a mutual distance d = 32 A calculated in the (a) local spin 
density approximation and (d) with a one-orbital mean-field 
Hubbard model, (b) and (e) panels show the same, but for 
the antiferromagnetic state. The resulting magnetoresistance 
in both approximations is shown in (c) and (f). 
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tion, here the electronic structure will be indirectly af- 
fected by the field even if its direct influence on the or- 
bital wave function is neglected. This change reflects in 
the spin-resolved conductance as shown in Fig. [Sja-b) 
for d = 32 A (the dashed hne represents the conduc- 
tance of the clean AGNR). The different total transmis- 
sion for the F and AF solutions, resulting from adding 
the two spin channels, results in a positive magnetoresis- 
tance (MR), MR — Gp ~ Gaf/Gf + Gaf^ at energies 
near the bottom and top of the conduction and valence 
bands, respectively [see Fig. EJc)]. Similar results are 
obtained (not shown) for different intermediate distances 
between H atoms. The right panels in Fig. [5] show the 
results obtained from the one-orbital mean-field Hubbard 
model for [/ — t — i eV. Apart from the recovery of the 
particle-hole symmetry, the results are remarkably simi- 
lar, validating the use of the latter model for transport 
calculations in hydrogenated graphene. 

One should note, however, that since the chemical po- 
tential lies in the middle of the gap, the energy ranges at 
which MR could manifest itself are not relevant in linear 
response transport for infinite AGNR's. A finite bias cal- 
culation may reveal the MR obtained at those energies, 
but this is a non-trivial task beyond the scope of this 
work. The application of a gate voltage is not a prac- 
tical alternative either since it implies a deviation from 
charge neutrality which would fill up or empty the local- 
ized states hosting the unpaired spins and kill the mag- 
netization. Instead, we propose to explore the possibility 
of MR at zero bias by considering finite AGNR's con- 
nected at the ends to conductive graphene. This is done 
in our calculations by considering a metallic AGNR with 
a narrower section in the middle of appropriate width 



FIG. 6: Tunneling magnetoresistance for four different atomic 
H configurations. Two H atoms in a head-to-head configura- 
tion located in the middle of an armchair graphene nanorib- 
bon of length L — 73.7 A with a width of A'^ = 9 (a) and 
A'' = 15 (b), and two H atoms located near one edge in a 
head-to-head (c) and in a tail-to-tail (d) configuration for a 
ribbon width of A'^ = 9 and length of L = 57.0 A. 



[see Fig. IHd)]. (Note that in the one-orbital mean-field 
Hubbard approximation AGNR's of width N — 3m — 1, 
where m is an integer, are metallic, being semiconduct- 
ing otherwise). In what follows and in the light of the 
previous results, we restrict ourselves to the one-orbital 
mean-field Hubbard model. 

In our proposed AGNR heterostructure the difference 
in the zero-bias tunnel conductance between the F and 
AF states is now responsible for the appearance of tun- 
neling MR (TMR), as shown in Fig. [6] Notice that 
unlike conventional TMR, where the magnetic elements 
are in the electrodes, in our proposal magnetism is in the 
barrier. Panels (a) and (b) in Fig. [5] correspond to H 
atoms placed in the middle of a semiconducting AGNR 
of length L = 73.7 A and width N ^ 9 a.nd N = 15, 
respectively. Both cases refer to head-to-head configura- 
tions. The obtained TMR changes sign with d, but it 
is always negligibly small. On the contrary, when placed 
near the edge [Fig. [BJc)] , the TMR is positive and reaches 
much larger values. This result is for an iV = 9, L = 57.0 
A AGNR. When the H atoms are now placed in a tail- 
to-tail configuration near the edge for the same AGNR, 
the TMR becomes negative and reaches values of up to 
100% [see Fig. [H^d)]. As a final example we show in Fig. 
[7] the TMR for the case where the H atoms are placed 
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FIG. 7: Tunneling magnetoresistance as a function of the lon- 
gitudinal distance between H atoms for three different arm- 
chair graphene nanoribbon widths. The If atoms are located 
on opposite edges of the ribbon. 

on the opposite edges of an L = 57.0 A AGNR for three 
different widths. Large (and negative) values are also 
obtained for the tail-to-tail configurations in the range 
of distances explored. As can be appreciated, the TMR 
logically decreases with the ribbon width. 

IV. PRACTICAL CONSIDERATIONS 

A critical assesment of the practical consequences of 
the results presented is due at this point: 

• First, the TMR results presented above have been 
obtained for a particular type of AGNR het- 
erostructure and rely on the existence of metal- 
lic AGNR's. Calculations where the relaxation of 
the atomic structure on the edges is considered re- 
veal that a gap always opens even for the nom- 
inally metallic armchair nanoribbons^, compro- 
mising the applicability of these nanoribbons as 
electrodes. Graphene-based metallic leads, how- 
ever, can be found in recent literature. For in- 
stance, zigzag graphene nanoribbons with passi- 
vated edges are metallioii. The physics described 
in this work does not rely on the edge termination 
of the nanoribbons and could be reproduced in fully 
passivated zigzag nanoribbons as well. Another 
alternative can be based on using partially un- 



zipped metallic carbon nanotubesai, as suggested 
in KeiM-. A third possibility consists in gating 
selectively a semiconducting AGNR as previously 
done for nanotubes^. 

• From the LSDA calculations we note that the en- 
ergy difference between the F and AF states, or 
exchange coupling energy, can be as large as tens 
of meV for short distances, particularly for head-to- 
head configurations. However, as far as magnetore- 
sistive properties is concerned, exchange coupling 
energies above 1 meV are of no practical use since 
the Zeeman energy gain per spin in a magnetic field 
B is OMlgnsB meV T^^ and fields higher than 10 
T are hardly accesible in the lab. The AF-F transi- 
tion is therefore practically forbidden below c? « 65 
A and d « 40 A in the examples shown in Fig. [Sja) 
and (b), respectively. At larger distances the fields 
necessary to induce the AF-F transition can be as 
small as needed, but also is the associated MR as 
shown in Fig. [6ja) and (c). This is not a problem, 
however, when the H atoms are coupled tail-to-tail 
as, e.g., in the cases where they are located on op- 
posite edges. There, as shown in Fig. [71 a TMR as 
high as « 50% can be obtained at reasonably small 
exchange coupling energies [see inset in Fig. [31(c)] . 

• The third caveat relates to the stability of the 
atomic configurations. Given the tendency for 
H atoms to form lowest-energy non-magnetic 
aggregates'^, the existence of magnetically active 
dilute ensembles of adsorbed H atoms requires cer- 
tain conditions. For instance, working at reason- 
able low temperatures prevents H atoms from di- 
fusing after adsorption^. Another possibility is one 
intrinsic to AGNR's: The binding energy is larger 
for H atoms close to the edge than in the bulk of the 
nanoribbon. Related to this is the fact that the mo- 
bility of H decreases significantly near the edge, re- 
ducing the possiblity of formation of lowest-energy 
non-magnetic aggregates near the edge^^. As we 
have shown, it is precisely the hydrogenation near 
the edge that favors the appearance of sizable MR 
and thus an actual experimental verification. 

V. SUMMARY 

In summary, we have shown, as a proof of princi- 
ple, that hydrogenated AGNR's at small concentrations 
can exhibit magnetoresistive properties without invok- 
ing purely edge-related physics. Our results, which are 
deeply rooted in Lieb's theorem, provide further evidence 
for the wide applicability of this theorem beyond the 
bipartite lattice Hubbard model for which was demon- 
strated. We have also shown that hydrogenation near 
the edge presents advantages with respect to bulk hydro- 
genation both from energetic and electronic standpoints. 
This aspect, especific to nanoribbons, might favor the use 
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of these ones for spintronics applications, as opposed to 
using large flakes of graphene. Although our conclusions 
are based on the simplest case of two H atoms, noth- 
ing seems to prevent magnetoresistance from occuring for 
emsembles of many H atoms. This, however, still needs 
to be supported by a more extensive statistical analysis 
which is beyond the scope of this work. 
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